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ABSTRACT 

Flowfield measurements of the blade tip vortices from a rotating 
helicopter rotor model were performed by three component laser- 
Doppler velocimetry (3D-LDV) and conventional (two component) 
particle image velocimetry (PIV). In general, the results are in good 
correspondence, but also illustrate the different properties of both 
techniques: LDV offers the capability of three-component measure- 
ments, whereas PIV captures the unsteadyness of the flowfield. 


1.0 INTRODUCTION 

With increasing use of civil helicopters the problem of noise emission 
of helicopters has become increasingly important within the last 
decades. Helicopter noise has been the subject of many research pro- 
jects<T. Blade vortex interactions (BVI) have been identified as a ma- 
jor source of impulsive noise. As BVI-noise is governed by the in- 
duced velocities of tip vortices, it depends on vortex strength and 
miss-distance, which itself depends on vortex location, orientation, 
and convection speed relative to the path of the advancing blade. 
Blade vortex interaction can occur at different locations inside the 
rotor plane depending on flight velocity and orientation of the blade 
tip path plane. Two types of BVI can be distinguished, (1) the angle 
between leading edge of the blade and vortex axis is about 0°, and 
(2) interactions at angles different from 0°. The noise emission of 
parallel blade vortex interaction is considerably larged 2 - 3 > compared to 
other cases of interaction. Investigations of the accoustic near and 
farfield 14 - 5 ) were based on the interaction of the blade with vortices, 
which were described by mathematical models. Information about 
the structure and strength of the real rotor tip vortices and their inter- 
action with the blade were not available. It is understood, that the study 
of these phenomena is of particular interest in progress towards quieter 
helicopters. 

In our investigations the vortical structures of the flowfield of a 
helicopter rotor model in a windtunnel has been studied by optical 
measurement techniques, as only non-in trusive techniques are capable 
of obtaining velocity data within the rotor plane. Measurements of 
local flow vectors at positions close to the rotor blade tips were 
performed by three component laser Doppler velocimetry (3D-LDV) 
and conventional (two component) particle image velocimetry (PIV). 
LDV measurements on helicopter rotor models have intensively 
been performed over the past two decades^ 6 ’ 8 ). The technique in 


general is described well in aerodynamic applications literature^ 911 ). 
The application of PIV to high speed air flows started approximately 
ten years ago( !2 > l3 >. Descriptions of multiple aerodynamic applica- 
tions are for example, given by Refs 14 and 15. PIV measurements 
on helicopter rotor models have been reported by Ref. 16. For the 
first time a comparison of LDV and PIV data acquired in a windtun- 
nel from the same rotor model has been performed. The blade posi- 
tion has to be chosen where vortex positions and their orthogonal 
interactions with the blade are well known and easily reproducible. 


2.0 ROTOR AND TEST MATRIX 

The LDV and PIV measurements have been performed on the 
helicopter rotor model of the Department of Aerospace Engineering 
(ILR) of the Technical University RWTH Aachen, developed and 
instestigated in detail by Ref. 17 (see also Fig. 1 and 2). The rotor model 
had a radius of 0-5 m and four NACA 0015 blades (chord length = 
54 mm) with rectangular tips. The rotor was driven by a toothed belt 
and an electric 65 kW engine and was fully articulated. It was 
installed in the open test section of the ILR Aachen low-speed 
windtunnel. 

The freestream velocity was set to U = 15-7 ms 1 and the rotor 
speed to/= 25 rev/s resulting in an advance ratio of p. ~ 0-2. The tip 
path plane was tilted by an angle of a TPP = -3° against the mean flow 
(forward flight), the collective part of the angle-of-attack was x> coU = 10°, 
and the cyclic pitch was x> cyc , = ±3-5° resulting in a trimmed condition. 
The Reynolds number, based on the blade chord length and blade tip 
velocity, was Re ~ 278 000. 

For the investigation of the blade tip vortex itself, 3D-LDV 
measurements have been performed at an azimuth angle of y ~ 0° in 
a plane outside the rotor (Fig. 3, position A). At this position the 
blade tip vortex, generated at y = 0°, is about 3 ms old and has rolled 
up completely so that there is no temporal change in the vortex structure 
itself. The blade load, and therefore the vortex strength at this azimutal 
position is similar to that at y = 90°. 

For the investigation of orthogonal blade vortex interaction LDV 
and PIV measurements have been taken at an azimuth angle of 
y ~ 90° (advancing blade side) in a plane perpendicular to the 
freestream velocity (Fig. 3, position B). 
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Figure 1 . 3D-LDV and helicopter rotor model in the open test section 
of the ILR Aachen low-speed windtunnel. 


3.0 3D-LDV SET-UP 

The 3D-LDV system of the DLR Institute for Fluid Mechanics was 
operated in a back scatter, off-axis models). A 6 Watt argon-ion 
laser was used as a light source of which the three most intensive 
laser lines (476-5 nm, 488 nm, 514-5 nm) were utilised to distinguish 
the different velocity components (Fig. 4). Each of the laser beams 
was divided into two individual beams with similar intensity, one of 
which was superimposed with a Bragg shift of 40 Mhz to reduce 
ambiguity. The beams were coupled into single mode glass fibres and 
were launched into the probe volume with individual transmitting 
optics. In order to provide a good resolution of all three velocity 
components, the angle between the optical axes was set to 30°. The 
effective size of the probe volume was approximately 0-25 mm in 
diameter and 1 mm in length. The tracer particles used were the 
same for both the LDV and PIV measurements: oil droplets, with an 
average diameter of less than 1 pm. As the receiving optics of the 



Figure 2. PIV set-up and helicopter rotor model in the open test 
section of the ILR Aachen low-speed windtunnel. 



Figure 3. Sketch of the paths of the tip vortices at an advance ratio of 
(p = 0-2, viewing direction parallel to the rotor axis, according to 
Beesten< 17 ) — not to scale. 


system operated in back-scatter mode, a sufficiently high data rate 
required the f-number of the receiving optics to be as small as possible, 
In this case a modified Schmidt-Cassegrain telescope with an 
aperture of 350 mm (f = 12) was used in order to gather sufficient 
scattered light from the particles. The received light, containing 
information of all three velocity components, was coupled into a 
multi-mode glass fibre, which transmited the light to a prism system 
in order to separate the different wavelengths. Then the light was 
converted into an electrical signal by three separate photo-multipliers. 

After digitisation of the signals, a fast Fourier analysis yielded the 
desired Doppler frequency. Only if the signals of the three components 
coincided, was the velocity data stored together with the actual time 
and measuring position. 

Based on the assumption of a periodic flowfield with respect to the 
rotor revolution, conditional sampling was used to get time dependent 
information of the flowfield. A trigger signal from the rotor axis was 
shifted with a time delay to any preselected azimuthal position of the 
rotor blade with an accuracy of 0-3% with respect to the rotor 
revolution. At this position, the data acquisition system clock was set 
to zero. Data acquisition was then started for a time window, corre- 
sponding to an azimuthal window or a time interval needed for the 
flow structure to pass through the probe volume. Usually data were 
acquired just within a small azimuthal window of 20° to 60°. 

Converting time information into spatial information can be done 
in different ways. In this case the transformation of the time informa- 
tion into an equivalent azimuthal angle was done by the following 
equation 



Figure 4. Sketch of the 3D-LDV system. 
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V - + 2n • / • t, where/is the rotor frequency . . • ( 1 ) 

The transformation of the time dependent information into a spatial 
information of the flow structure as for example the length I was 
done by analysing the convection speed of the flow structure 

l = U c ,„„-t ...&) 

Using the conditional sampling mode to acquire time resolved 
velocity data, which is later converted into spatial velocity data, 
leads to an inaccuracy due to the fact that flow structures may 
change during the time interval /, corresponding to the equivalent 
azimuthal angle. This has to be balanced against the significant 
reduction in measuring time. 

The LDV signals were typically analysed with 128 samples per 
burst, resulting in a frequency resolution of approximately 0-1%. 
Based on the performance data, and the geometry of the setup used, 
the measurement uncertainty has been estimated to be approximately 
0-3 ms 1 for the velocity components normal to the axis of the receiving 
optics and 0-75 ms-' for the component parallel to the axis. An 
additional error due to the velocity lag of the tracer particles has 
to be taken into account. This problem is discussed intensively in 
literature^ 9 ’ 20 >. However, in the case of the particles used for the LDV 
and PIV measurements — 1 pm or less in diameter — the velocity 
lag can be said to be sufficiently small with respect to other measure- 
ment uncertainties 121 ) to be neglected. 

A ‘position monitoring system’ gave access to blade motion parameters 
(lead-lag motion, pitching motion and angle of incidence) at a prese- 
lected radial position of the blade. 

4.0 PIV RECORDING SYSTEM 

During the last decade, PIV has increasingly been used to measure 
instantaneous flow velocity fields. This technique, in contrast to 
LDV, requires no conditional sampling. PIV allows the capture of 
the flow velocity in a two-dimensional plane within the flow within a 
few microseconds. It therefore allows for the entire velocity field to 
be obtained even in case of large cycle-to-cycle variations. The fact 
that the recording time necessary for the application of PIV (~ 12 ps) 
is small compared to the time required for one revolution cycle 
(40 000 ps), makes PIV an ideal tool for the investigation of the 
unsteady flowfields associated with rotor aerodynamics. 

During the past decade, a PIV system has been developed at DLR 
which can be operated under the rough environmental conditions in 
windtunnels< 22 >. This system utilises a dual-oscillator Nd:YAG pulse 
laser system with a maximum pulse energy of 2 x 70 mJ for illumi- 
nation of an area of up to 20 x 30 cm 2 of the flowfield. The recordings 
are taken with a 35 mm photographic camera and are analysed by a 
fully digital evaluation system. The cross-correlation software 123 ) 
was used to recover the displacement data. During acquisition, the 
laser power was limited to approximately 60 ml per pulse which was 
insufficient to saturate the sensor and produce adequately sized parti- 
cle images; many particle images thus only covered 1-2 pixels. As a 
result the measurement uncertainty was of the order of pm 01 pixel. 
The maximum displacements were limited to less than 3 pixels 
around the strongest vortex due to the short pulse delay of 12 ps. As 
a consequence, the measurement uncertainty was large compared to 
the measured displacements. Nevertheless, the velocity data uncer- 
tainty has been estimated to be below 0-5 ms- 1 for the measured 
in-plane components. As stated earlier, the tracer particles used were 
the same as for LDV. 

Using the synchronisation scheme shown in Fig. 5, the recording 
of the PIV images could be performed phase locked with the motion 
of the helicopter rotor model. Therefore, cycle-to-cycle variations of 
the flowfield at the same azimuthal angle of the reference blade 
could be investigated by analysing up to 100 recordings, each of 
which contained over 2000 independent velocity vectors. However, 
only the two in-plane components of the velocity vectors could be 



Figure 5. Synchronisation scheme of PIV data acquisition. 
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measured by means of the photographic PIV system. Examples of 
fully digital, non-photographic 2D-PIV measurements performed 
during the same measuring campaign are presented in Ref. 15. 


5.0 LDV RESULTS 

In Fig. 6 the tangential velocity profile of the tip vortex at \|/ = 0° is 
shown along a line through the vortex centre. For the measurement of 
the profile, the resolution was adapted to the extension of the vortex to 
be able to resolve the flowfield gradients. From this data the vortex 
core diameter has been estimated to be 4-8 mm, or 8-9 % of the blade 
chord respectivily, and a maximum tangential velocity of 18 ms- 1 . 

Figure 7 depicts the axial velocity profile of the tip vortex along 
the same line through the vortex centre as for the tangential velocity 
profile given in Fig. 6. For the axial velocity component, there is an 
increase in velocity of 4% to W ax /W lip = 1-04. 

As shown in Fig. 3, LDV measurements in a conditional sampling 
mode were collected in a plane perpendicular to the freestream velocity 
at the advancing blade at an azimuthal angle of \|/ = 90°. This yielded 
a three dimensional measurement grid, with the two coordinates 
obtained from traversing the LDV and the third coordinate derived 
from the transformation from time resolved velocity data to spatial 
velocity data. Associated with each of these grid points is a local 
three-dimensional velocity vector. By taking velocity data at each 
grid point for a time of 1-75 ms, corresponding to an angular window 



Y [mm] 

Figure 6. Tangential velocities along a line through the vortex centre 
obtained by 3D-LDV at \|/ = 0°. 



Y [mm] 


Figure 7. Axial velocities along a line through the vortex centre 
obtained by 3D-LDV at = 0°. 



X [mm] 


Figure 8. Velocity vector plot (2D) in a plane 0-2 chord length behind 
the blade tip as measured by LDV at y = 90°. The vorticity contours 
were obtained from the velocity data using finite differencing. 



time [ms] 

Figure 9. Temporal development of maximum vorticity of the blade tip 
vortex as measured by LDV at \\i = 90°. 


of A 1 ? = 10°, leads, with a resolution of 0-5°, to a number of 20 
measuring planes behind the rotor blade. 

First, vector plots were computed from velocity data to give a 
general impression of the flow velocity field. Figure 8 shows the vor- 
tex locations inside the measurement grid. At this azimuthal position 
a vortex can be observed, which has just been generated (age = 0 ms, 
position corrected by the tip displacement X = 7-5 mm, Y = 0 mm). 
The vortex generated by the 90° preceeding blade (age =10 ms, 
position corrected by the tip displacement X = 43 mm, Y = 5 mm) 
was located closer to the reference blade than expected (see also 
Fig. 10). Due to the location of the vortex generated by the advancing 
blade, an interaction with the blade occurs, which results in an 
additional vortex of opposite rotational direction (X = 25 mm, Y = —5 
mm). Since the time between the generation of the tip vortices and 
their measurement is small (max 10 ms) the cycle-to-cycle variations 
can be expected to be small. 

From each of the 20 planes, the maximum of vorticity of the newly 
generated blade tip vortex can be derived. This will give information 
on the temporal development. Due to conditional sampling, ten addi- 
tional planes in advance of the rotor blade are also available, which 
includes information on the maximum vorticity of the vortex generated 
by the 90° preceeding blade. Figure 9 shows the temporal develop- 
ment of the maximum vorticity of the blade tip vortex from ( = 0ms 







Raffjel et al Vortical flow structures at a helicopter rotor model measured by idv and piv 


225 


up to t = 12 ms. The vortex formation (I) can be described by a 
second order polynominal, while its maximum vorticity decreases 
inversely proportional to its age, i.e. 1 It (II). 


6.0 PIV RESULTS 

PIV measurements of orthogonal blade vortex interaction (position B) 
were taken at an azimuthal angle of \| / = 90°, again on the advancing 
blade side. The observation area was nearly parallel to the trailing 
edge of the blade and orthogonal to the axis of the vortices. Figures 
10(a) and (b) show two different instantaneous velocity vector fields 
from a set of 100 PIV recordings obtained at this angle. The origin 
was fixed to the tip of the trailing edge of the blade as it passed 
through the image plane with full aerodynamic load. Figure 10(c) is 
the velocity field obtained by averaging all 100 PIV data sets, and below 
(Fig. 10(d)) an estimate for the mean out-of-plane vorticity component 
is given. The tip vortex (A), which has just been generated (age ~ 0 ms), 
was located at X = 7-5 mm, Y = 2-5 mm. A tip vortex (B), previously 
generated by the 90° preceeding blade (age ~ 10 ms), is now located 
at X = 22 mm, Y= 25 mm. 

The cycle-to-cycle variations of these two vortices were small 
enough such that they are properly resolved in the average velocity 
field. Close inspection of Figs 10(a) and (b) does, however, show 
some variations in the shape of the vortices. Also, a third vortex 
structure (C), generated by the 180° advancing blade, can be 
observed. As this vortex intersects the rotor plane, at the same time 
the blade intersects the image plane, it is sliced by the blade (i.e. 
orthogonal BVI) such that only its remnants can be observed in the 
velocity field of Fig. 10. 

Due to small cycle-to-cycle variations, this third vortex is sliced 
differently each time, which results in the structure essentially being 
lost in the averaged velocity field (Fig. 10(c)). 

In Fig. 1 1 the tangential velocity profile of a tip vortex (age = 0 ms) 
from a single PIV velocity data set has been plotted along a line 
through the vortex centre. 

From these data, the vortex core diameter has been estimated to be 
74% of the blade chord at \\t ~ 90°, compared to 8-9% by the LDV 
method measured at \|/ ~ 0°. The maximum tangential velocities (±24 
ms- 1 ) are 30% higher than those obtained by LDV (±18 ms-'). Also, 
it can be seen from the velocity profile in Fig. 11, that the 
V-component of the measured velocity does not converge to zero at 
the outer edges of the observation area. Similar results were recently 
presented by Ref. 16. This effect, and the discrepancy between the 
velocity profiles shown in Figs 6 and 11, can be explained by the 
fact that the flowfield measured at \|/ ~ 90°, consists of different 
vortices and their interaction with the shear layer generated by the 
rotor blades. This made the determination of circulation of distinct 
vortices and a fit to analytical vortex models difficult. This is also the 
reason why further analysis of the vortices is based on peak vorticity 
instead of circulation. 

In order to investigate the uncertainty resulting from averaging the 
velocity field, the data of a number of instantaneous velocity vector 
fields (Fig. 10) have been averaged. This way, the amount of 
decreasing maximum tangential velocity and increasing vortex core 
size due to the averaging process can also be determined. 

Figure 12 shows the tangential velocity profile along a line 
through the vortex centre obtained by a spatial average over 36 PIV 
recordings. The maximum tangential velocity has been estimated to 
be ±19 ms- 1 and the vortex core size is 10-1% of the blade chord, 
both in good agreement with the LDV data. The maximum tangential 
velocity and the vortex core size are in the range of the results 
derived from the LDV data. 

In analogy to the available LDV data, the calculation of the temporal 
development of the maximum vorticity for the newly generated tip 
vortex can be performed. In this regard PIV recordings were obtained 
at different azimuthal positions of the rotor blade. From these the 
maximum vorticity in relation to the rotor azimuth was derived. 
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Figure 10. Two examples (a, b) of 100 instantaneous velocity vector 
fields as obtained by PIV at \|/ ~ 90°. Below the averaged velocity field 
(c) and (d) the corresponding vorticity estimate, oo z , is shown. 
(Contour levels are spaced at 400/s; dashed contours indicate 
negative values). Only portions of the 210 mm by 150 mm fields are 
shown for clarity. 
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As with the results derived from the LDV mesurements, the roll 
up process can be described by a second order polynominal, while 
the vortex dissipation is governed again by a 1 It dependence. 


7.0 DISCUSSION 

The 3D-LDV measurements yielded fundamental results concerning 
the structure of the blade tip vortices. In addition to geometric parameters 
like location of the vortex relative to the rotor plane and orientation 
of the vortex axis in space, the vortex core size, axial velocity, vortex 
strength, and vorticity distribution could be derived . Although blade 
motion and local velocity had been measured simultaneously, phase 
averaging and the rearrangement of the pointwise measured velocity 
data to form a complete velocity field causes some problems. Aperiodic 
flow phenomena, even of small amplitudes, can lead to the loss of 
information about the cycle-to-cycle variation resulting in an underes- 
timate of the actual vortex strength, and an overestimate of the vortex 
core size. 

However, recent LDV measurements in large aerodynamic facilities 
show that 3D-LDV data of high quality can be obtained with good 
spatial resolution^ 24 ). 

In spite of the difficult experimental conditions, high quality PIV 
data were obtained with sufficient spatial resolution. The time needed 
for data acquisition can be considerably decreased compared to the 
LDV measurements. A specially developed synchronisation scheme 
allowed the capture of instantaneous flowfield measurements at 
exactly the same phase angle of the rotor revolution and therefore 
the study of aperiodic features of the flow. Also with PIV, the location 
of the vortex relative to the rotor plane, vortex core radius, convection 
speed and vorticity could be measured. However, since conventional 
PIV measures only two velocity components, data, like the orientation 
of the vortex axis in space and axial velocity of a vortex, cannot be 
derived without changing the viewing direction. The use of a second 
camera in a stereoscopic arrangement would allow the measurement 
of all three velocity components without averaging data of different 
cycles. 

From the comparison between the LDV and PIV data, the differences 
in averaging and instantaneous data acquisition can be illustrated. By 
averaging a number of instantaneous PIV velocity vector maps, similar 
results to the LDA results are obtained, especially regarding the 
vortex core size and the maximum tangential velocity. The temporal 
development of vortex formation and its dissipation can be obtained 
from both the LDV and PIV results. 

Future investigations should be performed by instantaneous 
measurements of the three velocity components e.g. by stereoscopic 
PIV, which has already been used in aerodynamic researcW 25 - 26 \ By 
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Figure 1 1 . Tangential velocitiy along a line through the vortex centre 
as obtained by a single PIV recording at = 90°. 
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Figure 12. Tangential velocitiy along a line through the vortex centre 
as obtained by an average over 36 PIV recordings at = 90°. 
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Figure 13. Temporal development of maximum vorticity of the blade 
tip vortex as measured by PiV at \j/ = 90°. 


the measurement of multiple parallel planes, more information about 
the orientation of the vortices axes can be obtained. More detailed 
studies of tip vortices from stationary blades and their interaction 
with blades* 27 ' 28 >, should be performed in order to reduce the 
complexity of the vortical structures under investigation. 
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